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Abstract 



This study considers boiling binary mixtures of water with methanol or 
2-propanol at subatmospheric pressures. Liquid-phase equilibrium vapor 
pressures, binary phase equilibrium thermodynamic properties, heat trans- 
fer characteristics, and the critical heat flux (CHF) condition are determined 
for saturated pool boiling from a localized heat source while varying the con- 
centrations of methanol and 2-propanol in water. The heat source is an 
upward-facing copper surface submerged in a laterally-confined, finite pool. 
Low pressure boiling of aqueous mixtures provides a means of removing 
high heat fluxes while maintaining low surface temperatures. Small ad- 
ditions of alcohol to water increase the CHF condition above that of pure 
water. Higher concentrations of alcohol begin decreasing the CHF condition 
to that of the pure alcohol. While single-component correlations using mole 
weighted binary liquid thermodynamic properties have been shown to 
predict ideal binary mixture boiling behavior, they are unsuccessful in 
predicting the characteristics of aqueous mixtures. In this study, ther- 
modynamic properties of the non-ideal aqueous mixtures were either 
measured, predicted, or taken from the literature. The significance of these 
results to the use of binary coolants for electronics cooling applications will 
also be discussed. 
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1. Introduction 

The goal of practical boiling research is to predict boiling performance. Boiling performance 
consists of the superheat required for incipience, boiling surface temperature fluctuations, the 
fully-developed nucleate boiling heat transfer coefficient, and the maximum or critical heat flux 
condition. In applications where it is desirable to keep the temperature of a boiling surface low 
while removing the greatest amount of heat, reducing the saturation pressure may be a useful 
solution. However, reducing the saturation pressure may have an undesirable effect on other 
aspects of boiling. 

A reduction in the saturation pressure causes a corresponding decrease in the saturation or 
boiling temperature. Consequently, a superheat level is achieved with a lower surface tempera- 
ture. This approach is particularly useful when water, having a relatively high saturation tem- 
perature at atmospheric pressure, is used as the boiling liquid. Water is a desirable liquid since it 
has such a high heat of vaporization, high thermal conductivity, and is non-toxic and non- 
flammable. 

Boiling in sealed vessels is a typical application of subatmospheric pressure boiling. Heat 
pipes, thermosiphons, and some heat pump cycles may rely on subatmospheric pressures to 
provide low surface temperatures while moving significant quantities of heat. For example, it is 
often desirable to maintain a low temperature on the heated end of a heat pipe or thermosiphon in 
spot cooling of electronic components. Heat fluxes from current electronic components are ap- 
proaching 50 W/cm 2 . These fluxes are not easily handled by solid heat sinks. Phase-change 
heat sinks, which operate with a nearly isothermal interior, are becoming increasingly attractive. 
Low temperature operation of these heat sinks may be prescribed by creating a saturated liquid 
and vapor state in the vessel at subatmospheric pressures. Therefore the boiling occurs in the 
heated end of the vessel at a lower temperature. Knowledge of the boiling characteristics of the 
small heated surface is necessary to insure that steady and safe operating conditions are main- 
tained. 

The characteristics of pool boiling of water at low pressure are known to be much different 
from boiling at atmospheric pressure. Raben et al. (1965) investigated saturated nucleate pool 
boiling of water at subatmospheric pressures from a 3.81 cm diameter, horizontal heated surface, 
in an extensive pool. Cole and Shulman (1967) measured and correlated the effect of pressure 
on bubble departure diameters from a thin, 1.27 x 10.16 cm, horizontal zirconium ribbon. Van 
Stralen et al. (1975) studied nucleate boiling from an 11 cm diameter horizontal copper surface 
in an extensive liquid pool of water at subatmospheric pressures. In their experiments, bubble 
growth rates, frequencies, and departure diameters for different subatmospheric pressures were 
investigated. Single component boiling regimes of water and acetone at low pressure in a ther- 
mosiphon were examined by Niro and Beretta (1990). They used circumferentially heated 1.2 
and 3.0 cm I.D. by 20 cm vertical tubes. Van Stralen (1956) reported experimental CHF data for 
water as a function of pressure from a 0.02 cm diameter heated platinum wire. 
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Unfortunately, a result of boiling water at subatmospheric pressures is low-frequency bubble 
departure, creating surface temperature gradients and surface temperature oscillations. This 
phenomenon is seen in Figure 1 from the study of McGillis et al. (1991). Figure 1 shows the 
condenser vapor temperature versus time measurements for a heat flux of 5 W/cm 2 for water at 
three different subatmospheric pressures. Rises in the condenser vapor temperature indicate 
boiling activity. 

In some electronics cooling applications, it is necessary to keep the boiling behavior steady, 
with uniform heat removal. Pure water at subatmospheric pressures does not easily provide 
these attributes, particularly at low fluxes. Under some conditions, the boiling performance can 
be improved by adding a second lighter liquid to water. This would be classified as boiling of a 
binary mixture. Despite the numerous studies of binary mixture boiling processes over the past 
50 years, very little is known about the mechanisms responsible for their behavior. State-of-the- 
art review articles by Shock (1982) and Thome and Shock (1984) provide an extensive overview 
of prior multicomponent boiling research. There is also a dearth of information in the literature 
regarding boiling of water with organic liquids at subatmospheric pressures from horizontal sur- 
faces. Bonilla and Perry (1941) investigated combinations of water, ethanol, n-butanol, and 
acetone from a horizontal chromium plate at pressures ranging from 20 to 130 kPa. 

In this investigation, experiments were performed to determine if an aqueous mixture could be 
used as the working fluid in a two-phase thermosiphon used to cool electronics. This study 
considers binary mixtures of water with methanol or 2-propanol at subatmospheric pressure be- 
cause their properties are well documented, high purity levels are commercially available, and 
they are relatively safe fluids. Liquid-phase equilibrium vapor pressures, binary phase equi- 
librium thermodynamic properities, heat transfer characteristics, and the critical heat flux con- 
dition were determined for saturated pool boiling from a localized heat source. The concentra- 
tions of methanol and 2-propanol in water were varied. The apparent physical effects and sig- 
nificant findings are discussed. 
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Figure 1: Representative condenser vapor temperature-time 

responses in the thermosiphon from McGillis et al. [1991]. 
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2. Experimental Apparatus 

Figure 2 shows the experimental test section and system used in this investigation. The cop- 
per test section was machined to accommodate two cartridge heaters in the bottom end. The top 
half of the copper piece was milled to provide a long 1.27 x 1.27 cm square section. Within this 
section, 0.8 mm holes were drilled to the center to hold thermocouple wires. The copper and 
thermocouples were then cast in a low viscosity epoxy. The main body of the pool boiling con- 
tainer was made with 2.5 cm I.D. tubing. 
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Figure 2: Schematic diagram of test system and test section used in pool boiling experiments. 

In order to examine boiling at pressures below 10 kPa, various experimental designs were 
considered. Nucleate boiling is very dependent on cavity size, distribution, and wetting 
properties. The most difficult extraneous nucleation sites to control were at the interface be- 
tween the copper test section and the epoxy. This interface had to maintain a vacuum seal after 
thermal cycling. The Ablebond 342-13 epoxy that was eventually selected, adhered and sealed 
well to the copper. Once the epoxy cured, the top surface of the copper could be treated. Excess 
epoxy was milled down flush with the copper surface. The copper/epoxy surface was then 
finished with #320 emery paper and cleaned with alcohol. The epoxy surface was then bonded 
to the end of a glass tube to allow observations. The clear tube fit inside an O-ring fitting at the 
bottom of the condenser so that repeated assembly was simple. The condenser was made of a 
long section of copper tubing which had radial copper fins wound and soldered onto its O.D. 
Heat was removed from the fins with an air blower. 

The top end of the condenser tube was equipped with several sensors. A thermocouple probe 
extended down through the inside of the tube and could be positioned vertically to monitor fluid 
or vapor temperatures. A transducer measured the internal pressure of the thermosiphon. 
During initial start up of the system, a valve allowed a two-stage vane-type vacuum pump to pull 
the internal pressure down to very low values. The liquid was boiled during this step to remove 
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gasses from the internal volume. A second de-aerated fluid was added to the system via a septa 
and syringe to keep non-condensible gases out. Once the system was filled, the thermodynamics 
were verified. A computer and control algorithm varied the blower speed as needed, and proved 
able to keep the saturation temperature, T s , within 1°C during continuous boiling. This main- 
tained the internal pressure constant within 0.27 kPa. 

A datalogger recorded the temperatures and system pressure, with a sampling rate of up to 
1000 Hz. A linear fit of the measured temperature gradient in the copper section was used to 
calculate the test section heat flux and the heated surface temperature. Experiments and 2-D 
numerical analyses indicate that the heat losses from the test section were less than 6%, and that 
the average heat flux at the boiling surface, q", could be determined within 3%. Experimental 
uncertainties in the pressure and differential temperature measurements were 0.1 kPa and 0.1°C, 
respectively. 

Steady state for the entire experiment was determined by monitoring the temperature changes 
with time via the datalogger and computer. When a particular surface was extended to the criti- 
cal heat flux condition, the final, highest heat flux for which the system reached steady state for 
nucleate boiling was used as the CHF data point. 

Data presented in this paper were measured with a stationary liquid pool height of 7.1 cm. The 
pressures investigated ranged from about 7 to 83 kPa. Heat flux levels off the 1.27 x 1.27 cm 
surface were as high as 145 W/cm 2 . 
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The Boiling Process in the Thermosiphon 

The thermosiphon used in this study could operate in a variety of boiling regimes, from low 
heat flux, evaporative heat transfer, through nucleate boiling, to critical heat flux or dry out. 
When there was no heat supplied to the system, and the liquid and vapor were at thermodynamic 
equilibrium, the pressure in the system corresponded to the saturation pressure at the system 
temperature. The liquid became superheated with the application of very low heat fluxes. Heat 
was removed by evaporation at the liquid-vapor free surface and subsequent condensation oc- 
curred at the condenser walls. 

At higher heat flux levels, the liquid became superheated enough to allow bubble growth from 
cavities on the heated surface. As in any superheated system, the onset of nucleation in a cavity 
depends on: 1) residual gasses or vapor in the cavity; 2) the size, shape and material of the 
cavity; 3) the liquid and vapor thermophysical properties; and 4) the amount that the liquid is 
superheated. If all of the necessary criteria for bubble growth have been met, a bubble may grow 
and release from the heated surface. It is expected that when the bubble grows and departs from 
the surface, colder fluid will replace the highly superheated fluid. Since bubble growth depends 
on a sufficient superheating of the surrounding liquid, an appreciable time may be required to 
reheat the liquid in the vicinity of the wall to a superheated state and initiate subsequent bubble 
growth. This time is termed the waiting time, t w . 

The length of time from the beginning of bubble growth to bubble departure, t d , depends on 
how large the bubble must become for release to occur. This interval therefore depends on the 
rate at which the bubble grows to departure size. The departure bubble size is determined from 
the net effect of forces acting on the bubble as it grows on the surface. Surface tension holds the 
bubble down. For an upward-facing horizontal surface, buoyancy pulls the bubble up. If the 
bubble grows rapidly, the inertia associated with the induced liquid flow around the bubble may 
also tend to pull the bubble away. Bulk liquid motions may also produce lift forces on the bub- 
ble causing it to be pulled away. 

At higher heat flux levels t decreases, and the number of active nucleation sites and the fre- 
quency of bubble departure increases. As the flux increases, the boiling becomes more con- 
tinuous. 

The critical heat flux is not a boiling regime but the maximum heat flux attainable before the 
system makes a transition to film boiling. For most applications, a transition to film boiling is 
unacceptable because of the very large wall temperature excursions. 

3. Binary Phase Liquid- Vapor Equilibria 

The thermophysical properties of the fluid mixture must be determined before boiling perfor- 
mance can be predicted. For example, incipience, heat transfer coefficients and critical heat flux 
are heavily dependent on the fluid properties. State-of-the-art methods for predicting binary 
mixture boiling phenomena require interpretations of the mixture properties and the mixture 
phase equilibria characteristics. One's problems are not solved once the equilibrium binary mix- 
ture properties are determined. Non-equilibrium conditions and the transport of the properties 
have important effects. 
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The pressure and other thermophysical properties of a binary mixture can be determined from 
ideal and non-ideal mixture relationships. In both ideal and non-ideal mixtures, the chemical 
potential of component i in the liquid is equal to the chemical potential of component i in the 
vapor (see any thermodynamic textbook). The auxiliary function, fugacity, related to the chemi- 
cal potential, was created to simplify the abstract concepts of the chemical potential (see Praus- 
nitz et al. (1986)). It is shown that the fugacity of component i in the liquid is also equal to the 
fugacity of component i in the vapor. 

fi,v=A,l (1) 
The fugacity of component i in the vapor is related to the total pressure by 

fi, v = ^yi p (2) 

and the fugacity of component i in the liquid is related to the saturation pressure of pure com- 
ponent i by 

fi,l = Yi x i P pure,i (3) 
where (])• is the fugacity coefficient, y- is the activity coefficient, and y- and x i are the the vapor 
and liquid mole fractions of component i, respectively. For low pressure, (])•= 1 is a very good 
approximation. However, y- is not unity for non-ideal mixtures (see Prausnitz et al. (1986)). 
Substituting (2) and (3) into (1), it follows that 

i P = yi x i P pure,i ( 4 ) 
and using the fact that y i + yj = 1 and x i + Xj= 1, the saturation pressure of the mixture can be 
expressed as 

P mix~ x i Yi P pure,i + ^~ x i) Yj P pure,j @) 

Rearranging equation (4), the vapor pressure of the mixture is related to the vapor mole fraction 
by 

Y X P 
p _ ' i i pure,i 

mix ~ ^ W 

Experiments were performed in order to determine whether the non-ideal binary mixture be- 
havior agreed well with the current predictive methodology shown above. In these experiments, 
it was implicitly assumed that only small volumes of liquid were needed to occupy the vapor 
volume. Thus, the liquid composition was assumed not to change throughout the experiments. 
The saturation pressures for pure fluids and fluid mixtures were measured at different tempera- 
tures. Some measured pressures are plotted along with equation (5) in Figure 3. Figure 4 shows 
the liquid portion of the phase-equilibrium diagram for methanol-water at different temperatures. 
Measured values are plotted along with equation (5). It is quite apparent from Figures 3 and 4 
that the liquid phase equilibria can be predicted by the above models. 

Having determined and y-, the liquid-vapor equilibrium phase diagram can be determined. 
The phase diagrams for water/methanol and water/2-propanol are plotted on Figures 5 and 6 for 
a temperature of 38.9°C. Equations (5) and (6) provide the determination of the concentrations 
in the vapor and liquid at equilibrium for any system pressure and temperature. 
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Figure 3: Liquid saturation pressure versus liquid mole fraction for 
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Figure 4: Liquid phase equilibria for methanol-water. Comparison of 
experimental data with the non-ideal mixing model. 
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Figure 6: 2-propanol phase equilibria and surface tension versus 2-propanol mole fraction in water. 
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4. Discussion of Results 

Onset of Boiling or Boiling Incipience 

Bubble growth depends on surface characteristics, fluid properties, and available superheat. 
Bubble growth is expected to take place from an active cavity (cavity with residual gasses or 
vapor) when the surrounding liquid reaches the required superheat. The following estimate of 
superheat required for the onset of nucleation, AT 0NB , can be derived from static equilibrium. 

* T ONB = R Jp ldT)sat (7) 

Equation (7) indicates that the wall superheat required for the onset of nucleation is determined 
by the liquid surface tension, o (see Figures 5 and 6), the active surface cavity size, and the slope 
of the vapor pressure curve, dPIdT (see Figure 7). The value dPIdT for pure and binary mixtures 
can be calculated by taking the derivative of Equation (5). Fluids with a low surface tension, 
such as alcohols, will wet a given surface better and have a lower contact angle than higher 
surface tension fluids such as water. With a high surface tension fluid, the value of R c may be 
much greater than it would be with a lower surface tension fluid. As shown in Figures 5 and 6, 
the water/alcohol mixture surface tension is reduced from that of pure water. This change is 
extreme for the water/2-propanol mixture shown in Figure 6. 





0.2 0.4 0.6 0.8 1 
liquid mole fraction alcohol 

Figure 7: Variation in the slope of the liquid saturation curve, dP/dT, with mole fraction. 

As seen in Figure 8, for an average wall heat flux of 5 W/cm 2 and a pressure of 7 kPa, the 
boiling behavior for pure water is different from that of the x propano i=0.5 mixture. With the 
addition of alcohol to water, more and smaller bubbles are generated from the surface. 
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Figure 8: Photographs of boiling from a 1.27 cm square heated surface at 7 kPa. 

(top)Pure water. (bottom)x 2 _ propano[ = 0.5. 
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Boiling of Pure Fluids 

Figure 9 is a comparison of the nucleate boiling data obtained for water at two different sub- 
atmospheric pressures. A reduction in the pressure for a saturated water system shifted the boil- 
ing curve to higher wall superheat (T w —T s ) levels. However, the decrease in the saturation tem- 
perature associated with lower pressures more than compensated for this effect, resulting in 
lower wall temperatures for a given heat flux. The shift in the boiling curves may be attributed 
to a combination of effects. Lower pressures resulted in lower vapor densities and larger 
bubbles. The lower pressure increased the minimum superheat required for nucleation, resulting 
in a delayed onset. Also, greater fluid viscosities result from lower fluid temperatures, reducing 
the production of small eddies behind the growing and departing bubbles, thus reducing heat 
transfer. As can be seen in Figure 9, fluctuations in the wall superheat are more prevalent at 4 
kPa. 

The general increase of wall superheat with decreasing pressure is also seen in the pure 2- 
propanol and methanol boiling results shown in Figures 10 and 11, respectively. The pure 
methanol and 2-propanol boiling systems do not have discernible temperature fluctuations 
caused by intermittent boiling, however, both of these systems have a wall temperature over- 
shoot or excursion at heat fluxes less than 10 W/cm 2 . Both pure methanol and 2-propanol on 
copper surfaces are low contact angle fluids that wet the surface well. A well-wetting fluid can 
significantly reduce R c which, in turn, significantly increases AT 0NB , causing a large tempera- 
ture overshoot at incipience. 
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Figure 9: Effect of pressure on the nucleate boiling curves for pure water. 
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Figure 10: Effect of pressure on the nucleate boiling curves for pure methanol. 
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Figure 11: Effect of pressure on the nucleate boiling curves for pure 2-propanol. 
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As seen in Figure 12, for a given boiling system saturation temperature of 38.9°C, the pure 
water handled the highest heat flux but demonstrated unsteady boiling behavior at low heat 
fluxes. The wall temperature of the methanol system was the lowest and provided steady boiling 
behavior. The 2-propanol system operates with steady boiling behavior but has a relatively small 
critical heat flux and the highest wall temperature. 
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° pure methanol P=32.0kPa 
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Figure 12: Comparison of pure fluid boiling curves at the same saturation temperature. 
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Boiling of Aqueous Mixtures 

Figure 13 compares boiling curves of water and 2-propanol mixtures. The system temperature 
was maintained at 38.9°C. In each figure, the binary mixture boiling curve is plotted against the 
pure fluids at the same system temperature. Additions of water to 2-propanol decreased the wall 
superheat from that of pure 2-propanol. Additions of 2-propanol to water also reduced the fluc- 
tuations in wall temperature from that of pure water. This was due to the reduction in the bubble 
departure size. Although the average wall superheat did not significantly vary through the range 
of concentrations of 2-propanol in water, the magnitude of the critical heat flux did, reading a 
maximum in the vicinity of x 2—propanol = 

0.03. 

Figure 14 compares boiling curves of water and methanol mixtures. The system temperature 
was maintained at 38.9°C. In each figure, the binary mixture boiling curve is plotted against the 
pure fluids at the same system temperature. Additions of water to methanol reduced the wall 
superheat from that of pure methanol and reduced the fluctuations in wall temperature from that 
of pure water. The critical heat flux condition reached a a maximum in the vicinity of 

x 2-methanol = ®-^- 

Complex mechanisms govern the binary mixture boiling process. Binary mixture properties 
have multiple effects on the onset condition (Equation (7)) and therefore on the range of active 
cavities. The smaller the value of AT 0NB , the more sites of size R c will be active. An additional 
effect on the boiling process is the mass transfer resistance which takes place during the growth 
of the binary mixture vapor bubble. As a binary mixture vapor bubble grows from the surface, 
evaporation depletes the more volatile fluid (in this case alcohol) from the liquid-vapor interface. 
Because immediate diffusion or convection of the more volatile fluid to the interface is not pos- 
sible, the rate of evaporation is impeded. Consequently, the departing vapor bubbles are smaller 
and the amount of heat absorbed due to evaporation is reduced. 

By adding alcohol to water, the surface tension is reduced which increases the nucleation site 
density. The ONB condition is increased by a reduction in R c , causing an overshoot. The varia- 
tion of saturation pressure with temperature is increased slightly and should decrease the nuclea- 
tion site density. Mass transfer resistance reduces the bubble departure diameters and impedes 
the amount of energy absorbed, hence increasing the bubble departure frequency and reducing 
the fluctuations in wall temperature. 
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Figure 13: Boiling curves for varying concentrations of 2-propanol in water. 

Approximate fit to data for pure fluids plotted for comparison. T sat =38.9°C. 
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Figure 14: Boiling curves for varying concentrations of methanol in water. 

Approximate fit to data for pure fluids plotted for comparison. T sat =38.9°C. 
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Critical Heat Flux 

To avoid the risk of film boiling and accompanying high surface temperatures, the ability to 
predict the CHF condition is useful. The effect of pressure on the critical heat flux condition for 
pure fluids is shown in Figure 15. For the range of pressures considered, increasing the pressure 
increased the CHF condition. This is consistent with the CHF behavior observed for a variety of 
other pool boiling circumstances. A commonly used model which predicts the critical heat flux 
in saturated pool boiling for a surface of infinite extent is given by Zuber (1959). For low pres- 
sures, the Zuber critical heat flux is 

1 1 

^,z =0 - 131 Pv^/v^P/^)^ (8) 
Zuber' s analysis assumes the critical heat flux is attained when the large vapor jets leaving the 
surface become Helmholtz unstable. Notice that q" m z varies with V p v , so that with a given 
increase in pressure, there is a corresponding increase in the vapor density, and the critical heat 
flux is increased. An increase in the density of the vapor allows more energy removal per unit 
volume of departing vapor. 
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The critical heat flux condition of the aqueous mixtures considered in this study is not 
monotonic with concentration. Small additions of alcohol to water increased the CHF condition 
above that of pure water. Higher concentrations of alcohol began decreasing the CHF condition 
to that of the pure alcohol. Figure 16 shows CHF data for water/alcohol systems of this study. 
As seen in Figure 16, the data of this investigation could not be predicted by the basic Zuber 
model, particularly at low pressures. 
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Figure 16: Effect of concentration of methanol and 2-propanol in water on the critical heat flux condition. 
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5. Conclusions 

This experimental work investigated the boiling of aqueous mixtures from a small horizontal 
surface in a finite pool at subatmospheric pressures. The results support the following conclu- 
sions. 

1. In pure water systems, intermittent bubble departures cause large surface temperature oscil- 
lations at low fluxes. Large bubbles at low pressure cause significant downwash of cool liquid 
adjacent to the wall and subsequent wall temperature fluctuations. Binary mixtures reduce the 
size of the bubbles, thereby reducing the surface temperature oscillations. 

2. In pure alcohol systems, bubble sizes are reduced, relative to pure water, and temperature 
fluctuations are smaller. 

3. Water/alcohol mixtures at subatmospheric pressures provide a means of maintaining low 
wall temperatures without the large wall temperature fluctuations characteristic of pure water 
systems. 

4. Binary phase equilibrium properties are predictable and are an important element in under- 
standing binary mixture boiling. 

5. For pure water, methanol, and 2-propanol at subatmospheric pressures, the CHF condition 
for the small heated surfaces of this study deviates from the Zuber CHF predictions. The CHF 
condition for pure 2-propanol is small relative to pure methanol and pure water. 

6. For water/alcohol mixtures at subatmospheric pressures, the CHF condition for the small 
heated surfaces is not monotonic with concentration. The CHF condition is increased above 
those of pure fluids with small concentrations of the alcohol. 

The boiling behavior of water/alcohol binary liquid mixtures is not well understood or predict- 
able. Future research directions could include quantifying surface temperature fluctuations with 
fluid and pressure, correlating both nucleate boiling behavior and the CHF condition, and quan- 
titatively determining the actual vaporizing liquid concentrations. 
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